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PART I.- THE FLOW FIELD IN THE CORE REGION 


ABSTRACT: 

For the flow field of a free jet leaving a nozzle of circular cross 
section in a medium with straight uniform flow field, approximate formulas 
are presented for the calculation of the velocity distribution and the 
dimensions of the core region. The agreement with measured results is 
satisfactory. 
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I. INTRODUCTION AND SURVEY OF METHOD AND RESULTS 


Knowledge of the flow field of a free jet leaving a nozzle is of 
basic importance for practical application. 

Investigation of such a flow field is a problem of free turbulence. 

In theoretical research the following specialized cases of our 
problem have already been treated: 

(a) The mixing of two plane jets, the so-called plane jet boundary. 
These conditions are encountered in the immediate proximity of the nozzle. 

*"Zur Theorie der Ausbreitung eines aus einer Duse austretenden freien 
Strahls." Untersuchungen und Mitteilungen Nr. 8003, September 19^4. 
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(b) The spreading of a rotationally- symmetrical jet issuing from a 
point- shaped slot in a wall, the so-called rotationally- symmetrical jet 
spreading. This state defines the conditions at very large distance 
from the nozzle. 

In considering a free jet leaving a nozzle of circular cross sec- 
tion, we may subdivide the spreading procedure, according to an essential 
characteristic, into two different regions: 

(1) Region where a zone of undiminished velocity is still present 
(the so-called jet core). We shall call this range, which extends from 
the nozzle to the core end, the core region. For the immediate proximity 
of the nozzle the conditions of the plane jet boundary exist. 

(2) The region of transition adjoining the core region which is 
characterized by a constant decrease of the central velocity. This 
region opens into the region of the rotationally symmetrical jet 
spreading mentioned above. 

So far, there exists only an investigation concerning the core 
region (reference l) ; it is limited to the case where the surrounding 
medium is in a state of rest. 


Method and Results 

In the present paper, the spreading of a jet in the core region is 
treated for the general case where the surrounding medium has a straight 
uniform flow field (or, respectively, where the nozzle from which the 
jet issues moves. at a certain velocity through the surrounding medium at 
rest) . 

The theoretical investigation is based (reference 2) on the more 
recent Prandtl expression for the momentum transport 


e(x) 


Kb(x) 


^’max ” ^in 1 


One then obtains in the rotationally symmetrical case the following 
equations: . 


d(ru) d(rv) n 

+ “sr - = 0 


Continuity: 
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Momentum transport: 


— chi — du 

u + v 5? 


= € 


-J1 + I 

3r 2 r 


5? 


where 


e = KbCxX^ - u 0 ) 

Uq_ = velocity at the jet core 

Uq = velocity of the surrounding medium 

( u l> «o) 

With reference to the present problem 



we introduce, instead of r, 

as independent variable. 
Continuity 



We obtain: 
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Momentum transport 


u(x 


du 



— du €(x)/s 2 u 



Velocity distribution in the Core Region 


We limit our considerations at first to small disturbances of the 
flow field; that is, to relatively small differences in velocity 

U ° small quantity The partial differential equation for the 
momentum transport may then he linearized 


( U 1 - * 

\ U;L 


du , \/d u , „ S 2 u 


where 


c(x) = Kb(x) 


u i - U 0 

ui 


(it should be noted that by the transformation |i = 
is transformed into the equation 




this equation 


,d 2 u <3u 1 du 
£^2 + Sil n " ^ 


= 0 


which represents a heat conduction equation.) With q 
of r one obtains from the equation of momentum 


r - Tr 


instead 


c) 2 u 


St) 


- , du 
2 + 5) 


,♦2 

X 


+ 3 


^y 


du x 2 

SE v € ( xT 


= 0, 


n + ? / 0 
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This equation is a linear partial differential equation of second order 

of parabolic type. For the plane case ^ — y mj this results in the 
equation 


d 2 u 

dn 2 


+ 1^/* \ =0 

dn 


with 


n . b(x) 

= hm — - — — 


therein c = lim 

— =0 

r o 


i x 
e(x) 

r, 


0 


- -(^) 


and is to be regarded as a function of 


U 1 - ^0 
U 1 


With the boundary conditions taken into consideration, the integration 
yields 


with 


0 O = 


1WH 5 ) 


We now obtain an approximate solution of our ^problem by generalizing 
the plane velocity distribution and setting up the following formulation: 

z-M*- 1 ) e ‘ El "° 2]2 + °J + K 1+ W 
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With 



we obtain a function which corresponds to the exact solution for small 

— as well as for large positive r\, thus in boundary zones of the 

r 0 

region of integration as well as in the interior of the region along the 
jet t\ = 0. 

If we now consider larger disturbances, the solutions obtained for 
small disturbances are to be regarded as a first approximation. 

U 1 " U 0 

For the plane case the solution for arbitrary already 

1 

exists, compare Gortler (reference 3). It is found that, purely with 
respect to shape, even the first approximation represents^ very good 
approximation. The velocity distribution calculated by Gortler still 
shows an uncertainty insofar as u(ti + a), with u(t|), also represents 
a solution. This uncertainty here may be eliminated, because for the 
jet core vanishing of the transverse component v is required. There- 
with the initial profile of the velocity distribution for arbitrary 

Uf - Uq is then unequi vocally fixed. 

U 1 

If we limit ourselves, with respect to shape, to the first approxi- 
mation, the initial profile is 



where 
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and 


a 0 


/ U 1 - ^ 


! 2kc [ 


Un 


For the further development of the profiles starting from this initial 
profile the regularity found for small disturbances is then taken as a 
basis 


where 


JL = _L 

U1 \]n \U1 








+ 02 


- Uq 

For - — >0, this function is transformed into the approxi- 

ui 

mation function constructed for small disturbances. How far it may be 
considered an approximation in the region for arbitrary disturbance is 
not investigated in more detail. 


The functions appearing in the integral CTp(x), C2(x) result from 
the approximation calculation for the dimensions of the core region, 
carried out on the basis of the momentum theorem. 


Calculation of the Transverse Component 

The transverse component v of the flow is determined from the 
continuity equation 


v " ' rf ( r s) dr < r * °> 


or v 


1 





+ 


12 

X 



du 


~ T) 


du\ 

5F/ 


IdT) 


approximate function being substituted for u. 


respectively, with our 
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The integration constant is determined from the requirement that 
at the jet core the transverse flow component vanisjhes. 

In order to avoid complication of the calculation, rectilinear 
course of the mixing width b(x) is assumed. This assumption proves 
approximately correct as results from the calculation of the dimensions 
of the core region. 


Dimensions of the Core Region 

The dimensions of the core region (jet core and width of the 
mixing zone) are calculated according to a formulation of the momentum 
theorem 



(= rTxy) 




du 

57 


indicated by Tollmien (reference 4 ) . 


The occurring integrals as well as the ^ defining the shearing 

stress are determined approximately with the course of the velocity 
distribution assumed rectilinear 



Then there result for the limiting curve d(x) of the jet core and the 
width b(x) of the mixing zone two ordinary differential equations of 
the first order which can be reduced to one equation 


I - 


f(x) dx 


This integral can be represented with the aid of elementary functions; 
however, for simplicity its calculation here is performed by graphical 
method. 


k appears as the only empirical constant which results by compari 
son with measurements given by Tollmien (reference 4) as k = 0. 01576 . 
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Comparison with Measurements 

In order to carry through a comparison between theory and experi- 

up - Uq 

ment, a measurement for the case = 0.5 was performed with the 

test arrangement described in reference 5* 

The comparison with the theory offers satisfactory results if one 
takes into consideration that the effective radius of the nozzle flow 
referring to a rectangular velocity distribution is different from the 
geometrical radius. 


II . CALCULATION OF THE FLOW FIELD 
(a) Velocity Distribution in the Core Region 


We base the theoretical investigation on the more recent Prandtl 
expression for the turbulent momentum transport 


e(x) = Kb(x) | ^ - u^l 

where K = dimensionless proportionality factor, b - measure for width 
of the mixing zone, and u = temporal mean value of the velocity. 

We have at our disposal, for calculation of the flow field, the 
continuity equation and the momentum equation for the main direction of 
motion, which read in rotationally- symmetrical rotation 

Continuity 


d(ru) , d(rv) ^ 

“5T~ + ~5r = 0 


Momentum transport 


— chi — c)u 

U c>x Sr 


e(x) 


d^u + 1 du \ 

Sr 2 r * T ) 


where 


e(x) = ttb(x) ^up - uq) 
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u.]_ = velocity of the issuing jet 

Uq = straight uni form velocity of the surrounding medium 

> UQ 

We may integrate the continuity equation by introduction of a flow 
potential t 


- St 
m = 57 


- St 
^ = ' 57 


The momentum equation then is transformed into 

St S 2 t 5td 2 t 1 St St _ / J 5 3 t <£t -l|t\ 

57 57~57 “ 57 ^2 + r 57 57 " €W \ r Sr 2 ‘ Sr 2 r 57 ) 


where e(x) = Kb(x) 




- "o' 


U 1 


if we make the velocity dimensionless by 


division by up. According to a method applied by Gortler (reference 3) 
we set up for t the expression 

,2 


t = to + 


U 1 - " 0 

U 1 


fVL-i ~ 


tn + 


"0 


U1 


t- 


2 + 


'"1 - "0 

ui 


developing t in powers of the parameter 
the potential of an undisturbed flow ^u^ = u^j; thus 


Therein tg is 


8 *o 

ST ' ra i 


sr 


= 0 
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If we enter with this formulation into the differential equation, we 
obtain 




If one arranges according to powers of 
of differential equations for \|r^, 


/ U 1 ~ up 

\ up 


) 


one obtains a series 


For ^ 


^1 (A) ') + A (**o\ . A (**o\ 

dr \dx dr/ dx dr \dr / ^2- \dx ) 


/A\ 1 ^1 ( h %\ 1 ^1 Mo \ 

dx \^ r 2 / + r dr \dx ) + r dx \dr J 


= e(x)| 


d 3 *! 

dr^ 


d 2 ^ 1 dt x \ 

d^" + r 'ST j 
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or, taking 


»0 

* ru i 


into consideration 


^0 

sr = 


0 



dx Sr 


r 



a^i x stA 

i^ + ?5r / 


etc . 


On the Theory of Small Disturbances 


In the following, we shall limit ourselves at first to small 
disturbances of the flow field; that is, relatively small differences 

/Ui - Uq \ 

in velocity ( small quantity J, 

The velocity field is then defined by the flow potential ^1 • 
Since = ru, the above equation for may be written as 

follows : 




+ r 


djl\ 
hr 2 1 


e (x) = *b(x)( X Ui — ) 


Therewith we have attained for small disturbances a linearization of 
the equation of motion. 

(it should be noted at this point that by the transformation 
^ = r ■ our equation is transformed into 
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With reference to Reichardt's discussions (reference 6), it is of 
interest to point out that this equation is of the type of a heat 
conduction equation.) 

In view of the conditions existing in our problem 



(tq = nozzle radius, x = distance from the nozzle in direction of the 
jet axis) , we introduce instead of r the variable q = 


r - Tr 


coordinate transformation yields 


du _ du Sq _ 5u 1 
Sr Sq Sr Sq x 


This 


/ r=const 
thus the equation 


/du \ _ /Su\ 

\5£yL_ + '^*'q=const 


<3u Sn du du q 
Sq Sx ” Sx ” Sq x 


* r o>(s - - e(x> 


du 1 / . 1 

* x + r °) 5? ? 


— i 4 0 


or, respectively, for q + 


x 
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<3j3 

*n 2 


+ 



1 



x 


+ 11 ikj 


Su x 2 
3 x € (x) 




This equation is a linear partial differential equation of the second 
order of parabolic type. * 

The solution of this differential equation is fixed unequivocally 
by the initial condition that for ^ — > 0 the velocity distribution 
of the plane jet boundary appears. 

We first derive (for small disturbances) the velocity distribution 
of the plane jet rim. 


For — ->0 we obtain with the expression u(q) the equation 

r 0 



with 


x 

e(x) 


r 0 U 1 " "0 

— -=u KC ; 

r 0 Up 


b(x) 

X 

With the boundary conditions 


Therein c ~ lim * — ~ and is to be regarded as a function of 

-^=0 
r 0 


U 1 - ^0 
U 1 


u- 


up for T) } - 00 

Uq for q > + 00 
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taken into consideration, the integration yields 



Turning now to our problem, we can expect great difficulties in con- 
structing the exact solution. We limit ourselves therefore to forming 
an approximate solution. For this purpose we generalize the plane 
velocity distribution (the initial profile) and set up the following 
expression 


JL i_/Ho 

“i " \m W 


r 


[j f i(x)^+a2( x )l -[ai(x)T|+a2( x 71 


d|^CT 1 ( x) T| + a 2 (xj] 


1 

2 



This formulation insures at the outset a reasonable shape of the 
approximation solution. 

For <* 1,02 there immediately result, because of the initial 
condition, the requirements 


lim 




lim crgCx) = 0 

£ ^0 


Now the following equation is valid: 
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Accordingly, we put 


<Jl(x) 



g(x) 


(*l 

-r 


u o\ b 


U 1 


x 2cn 2 


x 

Furthermore we take care that our approximation statement for small — 
yields the exact solution. This will he the case when the ^j x ^ 


of the approximation statement agrees with the 

x 

from the differential equation for — = 0. 
According to the differential equation: 


(!) 


r o 

to be calculated 


du 


r 


du 

1 1 n —T—r " 

L 2 + ^ 

r 0 ^ e(x) 

T| + 

1 X 

/ 


Thus 


'3u \ 

&): 


-=0 


r 0 


= lim 
r 0 


_x 

r 0> 


= lim 




r 0 


=0 


— >0 


e(x) 

d 2 ^ , du 

-X 

1 + n * . 

2 < 
x c 


r 0 ^ e(x) 

T| + 

1 X _J 

— 

e(x) 

X 

1 

2oi 2 


1 

Su 


r 

0 2a i c 

! l\p 

L \ r ° 

> ® 


or, with 
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We now enter into this equation with our approximation expression; that 
is, we put (except for a common factor) 



We then consider the relations 



The last relation signifies that the width b of the mixing zone is, 
in the proximity of the nozzle, of rectilinear character, an assumption 
which seems justified considering the fact that we approach, in the 
proximity of the nozzle, the conditions of the plane jet boundary. 

We then obtain for the left side of the equation 
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for the right side 


, 8u 2„ 2 

Z2 Sn n 1 


^->o 2a i l + Ur 

r o 1— 


Equating yields the equation 




D) e 


( OqT] ) ‘ 


or respectively. 


ag'Co) = 2(^0 ” ^ a 0 0 2 * ^ ^ l| 


ff 2 ’(0) = 


This results, in a 2 


■ 


for small — . 


This guarantees first of all that our approximation expression 

for — > 0 represents the exact solution. 

r 0 . 

If we enter with the approximation expression thus constructed into 
the differential equation, we recognize immediately that the latter (due 

4-„ 4-v.„ ^ 2 ] j_ s satisfied also for T] * m 


to the factor e 
arbitrary — \. 


00 (and 


Thus our approximate expression with a^, <J 2 fixed in the above 

manner yields a funct ion which corresponds in boundary zones of the 
region to the exact solution. 

As to the behavior of our function in the interior of the region, 
it is found that the function in case of suitable "continuation" into 
the interior of the region satisfies the differential equation along q 
• ’ / 
For t| = 0 the differential equation reads 
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If one enters with the approximation expression and considers 

g(x) _ _1_ 

X 2 a l 2 

one obtains 



As solution one obtains 


x 



X 

For small — one has again 
r 0 



Therewith we have obtained for small disturbances the following 
approximation function 
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where 



To sum up: This function satisfies the differential equation with 

the initial conditions prescribed for small — as well as for large 

r 0 

positive qj in the interior of the region it satisfies the differential 
equation along the jet q = 0. Therewith we have constructed an approxi- 
mate function which in boundary zones of the region of integration and 
in its interior along the jet q = 0 is to be regarded as exact solution. 



Up - uq 

Let us now consider larger disturbances not a small 

u l 

quantity . 

First, we shall treat the problem of the initial profile. 

Gortler' s calculations (reference 3) showed that even the first 
( Uf - Uq\ 

approximation 1 for small 1 represents, purely with respect to 

shape, a very good approximation. This applies, however, only to the 
shape of the distribution curve - not to its position. The velocity 
distribution calculated by Gortler is unequivocally fixed by the 

__ Ui - Uq „ 

arbitrary requirement that u(o) = g • However, Gortler 

points out that with u(q), v(q), the equations u* = u(q + a), 

v* = v(q + a) - au(q + a) also represents a system of solution. But 
this remaining uncertainty is here eliminated by the fact that for the 
jet core the transverse component v must vanish as follows from the 
continuity. 

If u(q) , v(q ) is the velocity distribution calculated by Gortler 

U X - Uq 

which is characterized by u(0) = g , the quantity a must there- 

fore be determined in such a manner that v^ - au-^ = 0 which 
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V 1 


a = — . Taking Gortler's calculations as a basis, one obtains 
U 1 


yields 
in first approximation 



Therewith the initial profile for all — — is unequivocally 

determined. If we base the shape representation on the first approxi- 
mation, the initial profile is 


where 



and 


Hh*) 

For the further development of the profiles in the core region, 
starting from this initial profile, we take as a basis the regularity 
found for small disturbances. 
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where 

Ui - Uq 

t}* = cr-jTi - O.36 — + a 2 

U 1 

with the terms a^x), a 2 (x) determined before. 

This function therefore yields the initial profile in first approxi- 
mation. How far it may be regarded as approximation in the region is 

Ui - Uq 

not investigated in more detail here. For >0 it is trans- 

it 

formed into the approximate function found for small disturbances. 

Our approximate function generalized to arbitrary disturbances 
therefore reads 


where 


with 




NACA TM 1311 


23 


where 


x 

a , , r a ia) 

°1 ^ a 

Thus the curves result from one another by similarity transformations. 

Calculation of the curves requires, furthermore, knowledge of the 
functions o^x), a 2 (x) and, respectively, of the mixing width b(x) 

and the constant k. These quantities result from the approximate 
calculation (carried out with the aid of the momentum theorem) for the 
dimensions of the core region. 

the parameter values 

= 1.0} 0.8; 0.6; 0.4; 0.2 

— calculated for — =0 and the core end. 
U 1 ’ ■ r 0 

In figure 2 the functions a^x) and cr 2 (x) are plotted for the 

x 

parameter values named above, as functions of — up to the core end. 

Calculation of the Transverse Component 

The transverse component v of the flow is determined from the 
continuity equation 

8(ru) 8(rv) 

“§x + ST~ = ° 

and, respectively 



Figure 1 contains for 

U 1 - u 0 
Ul 

the velocity distributions 
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The integration constant is determined from the requirement that in 
the jet core the transverse component v must vanish. 

As the lover limit ve choose accordingly the determined by the 

bounding of the jet core (concerning the dimensions of the core region, 
compare next paragraph) . 

In order to avoid complicating the calculation, a rectilinear course 
of the mixing width b(x) was assumed. This assumption is approximately 
correct. (Compare fig. 11.) 


For the velocity distribution — ve substitute our approximate 

U 1 

function. The performance of the calculation (appendix no. l) yields 
the following final formula. 



S 
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+ a 



+ 2 Or 


0.36 










0 


)} 


Therein 

□ - (’ll - 0-36 ^^2 * o 2 ) 

n 0 -(<> 1 '%- °-36 Ul ^ * 0 g ) 



and F l'f] Snd F l’-[] 0 ’ respectively, signify the values of the error 
integral taken at the points [] and -CJ 0 > respectively. 

In figures 3 to 7 the distributions of the transverse component for 

a section (~~) =0.1 near the nozzle and a section of 3 /k of the core 

length are plotted for the parameter values — _ = l.o, 0.8, 0.6, 0.U, 

0 . 2 . U 1 

» 

U 1 - Uq 

In the case = 1.0 there are shown, moreover, the distri- 

butions for the sections l/k of the core length and the core end itself.* 

(Remark: The transverse components calculated for the core end seem 

to yield too small values of the approach flow; the reason is that the 
poor approximation of the velocity distribution, an essential charac- 
teristic of the Prandtl expression^ in the boundary zones takes the more 
effect in the calculation of the v component the more one approaches 
the core end.) 

Ul - Uq 

For small — — , the transverse component becomes very small 

(note the different scales in the various representations) . 
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(Id) The Dimensions of the Core Region 



The dimensions of the core region are defined by the limiting curve 
of the jet core d(x) and the width of the mixing zone b(x) or, 
respectively, the outer limiting curve of the latter b(x) + d(x) . 

According to Kuethe's procedure (reference l) we take as a basis 
the theorem of momentum in Tollmien's formulation (reference 4). 



X 
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If one marks off a control area in the indicated manner, one 
obtains in the rotationally symmetrical case 



Uq = velocity of the medium surrounding the jet. 

According to the more recent Prandtl expression 


T xy ~ Kb ( x ) ( u l - Uq) I? 

Thus we obtain, if we, furthermore, take the limits of the mixing 
zone into consideration 



p d+b 

(u - UqJ vr - ^ 

J u(u - Uq^t dr 


" r K b( Ul - u 0 ) ^ 


According to the existing conditions we transform (according to Kuethe) 
with 


Then 


*1 


r - d(x) 
b(x) 


- d (x) 
b(x 


r = bT) + . d 



If we make, in addition, the assumption that u depends only on t) 
not on x, there follows 
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For v we finally insert the continuity equation 




dr 


or 


XT 


v = 


(nb + d) 


^ - Ti (l) 2 bb’ - Ti(bd* + b'd) - dd' 


dq 


For approximate calculation, we write for the velocity distribution the 
sample expression 


u - uo 

1 

r-i 


u i 

l A 

1 - ') 


V 


r - d(x) 
b(x) 


This expression, which may be regarded as a first rough approximation- 
for the velocity distribution, will probably lead to not too large errors 

' for the integral calculation. The value ^ determining the shearing 

stress also will probably result in a usable approximation for the 
central region of the mixing zone. 


The result is 


8^ 

u l 

sr = ■ 


1 U 1 - "0 
b up 


_d_ 

u 

- Uq^ 

dq 


V Up / 


fUi - 


V 


U 1 


(1 - tl) 


^o foi - "o \ 

u l\ U 1 ) 


(l) We now put r = 0. The momentum theorem is then transformed by 

integration into the form of the theorem of conservation of momentum 



const . 
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or 


nd(x) rd+b _ . 

J u l( u l - ^i 1 * dr + J w(u - UQjr dr = u-^u-l - u 0 ) 


r 0 

2 


,2 0 1 


u l( u l - u o) u ( u - u o)( 7 i b ' + d ) b dT l = u l( u l “ ^o) — 


one obtains 


1/-2 2 \ 
2 \ - r 0 ) 


2\ 

r 0 ) 

♦»7 

J o 

u_/u 

ulv 

- u 0\ 
U 1 J 

u 

/Ui - Uq 

7- 

q) a 

U 1 

\ U 1 

+ b 2 

7ui - UqN 

u + 

% 1 


\ U 1 / 

/ 12 

u^ T> 


-1 _ 


carries out the integration, 


^ J / U l ' "0 1 "0 l\ 

( U 1 3 + u x 2) " 


0 


or 




+ bd 


1 _ lf u l - "0 


3\ u x 


,2 2 

+ d = r 0 


(2) In order to obtain a second equation between b and d, we 
put r = Tq. _ 

If one performs the somewhat .lengthy elementary calculation, one 
obtains finally (compare appendix no. 2) 


bb'< 


/* 

r~ 

1, — 

-\ 

U 1 - uo 

- U r ° ' 

i( + 1 

lFO-'jf 1 

ui 

k 

Z\ b 

s' 

) + ZJ 

,3V • b / " 3 


< Vd + + |J + i 


- d 


If* 


*1 ~ **0 
ui 
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The theorem of conservation of momentum reads in differentiated 
form. (Compare (l).) 



(h'd + bd') 


1 l( u l - M 

2 * Ul ~) 


+ dd' 



0 


By addition of the two equations one obtains 



(3) We now proceed to determine b and d from the two equations 
obtained. We replace b in the second equation by the expression for 
the function which we obtain by solving the first equation with respect 
to b. 
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Substitution then yields 



The evaluation of the integral could, in itself, be carried out 

by analytical method since the integrand is built rationally in — 

r 0 


and 
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a square root. However, the breaking up into partial fraction which 
has to be done in this procedure is very troublesome. Hence it is 
advisable to perform the evaluation graphically. 


For — =1 the integrand 

r 0 


h -Pir) 


fr 


assumes the indefi- 


nite expression The limiting value is 


lim f 3 

- d -_*l ' 


"i • u o 


lP> + =5 




2 / ~2 


t \ C 

I— I 

\ a 0 / 


/ A ^2 

i p. - u oj \ a ° a 0, 


£)' < 


Zk 

a 0 / 


If — was determined, analytically or graphically, as a function 

of — , b(x) results from 
r 0 


b d 

— = a 0 — + jn - a 2 




The relation 


\dxJjc__Q 


u i ^ 

Ut /V u arJ 


r 0 


lim 

d 1 
>1 

r 0 


f l * 


Ui - Uq 


Ui 


(which by comparison with measurements on the plane jet boundary may 
serve for the determination of k) also is of interest. 


The symbol k appears as the only empirical constant. 

With the measured results on the plane jet boundary with zero outer 
velocity (given by Tollmien (reference 4)) as a basis, there results 
with 



0.255 
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0.255 = - k(-3) 


0.185^ 


K = 0.01576 

Example a : 

In figure 8 the dimensions of the corresponding core region are 
represented for the parameter values 


ui ~ up 

Ui 


= 1.0, 0.8, 0.6, 0.4, 0.2 


x k 

Figure 9 contains the core lengths — , figure 10 the mixing 
b k r 0 u L - Uq 

widths — at the core end as functions of . 

r 0 U 1 


Figure 11 shows the mixing widths — for the various parameter 

U 1 - Uo x r ° 

values of as functions of — . 

ui r 0 

Figure 12 represents the angle of spread of the respective mixing 
region c - (g) ■ 

r 0 


Figure 13 represents ctq = 




as a function of 


1 ~ ^0 
U 1 


U 1 - Uq . Up - Uq 

with Tollmien s value c = 0.255 for = 1 being the 


£ 


U 1 

”0 


U 1 

defining quantity. 

Figure 14 shows for the medium at rest 

/db\ ' 'l 

as a function of c = K— ) . Figure 15 shows o 0 = — — 

\<ix/x_ = o \J2kc 

/db\ r ° 
function of c = |-r-| 

W4, o 


= 11 the quantity K 


as a 
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Figure 1 6 , finally, contains the limiting value 



necessary for calculation of the integrand in 


d 



III. COMPARISON WITH MEASUREMENTS 


Measurements on a free jet issuing from a nozzle and spreading in 
moving air of the same temperature do not exist so far. 

In order to test the theory by experiment, a measurement for the 

Ul - Ug 

case = 0-5 was performed at the Focke-Wulf plant. , 

U 1 

The measurements were carried out with the test arrangement with 
the 5 millimeter nozzle described in reference 5. A certain experi- 
mental difficulty was experienced in producing temperature equality in 
the two jets; it was achieved by regulation of the combustion chamber 
temperature with the test chamber pressure p^ and the probe pres- 
sure p kept constant. However a perfect agreement of the jet tem- 
peratures could not be accomplished inasmuch as the temperature measure- 
ment performed with a thermoelement is rather inaccurate in this low 
region. 

The test data were: 

Outer jet: 

Static pressure = -100 mm Hg 

(Measured relative to atmospheric pressure) 

o 

Room temperature tg = 20 

P 0 = 75^.5 mm Hg 


Barometer reading 
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Inner jet: 

Total pressure p a = 3^0 mm. Hg 

(Measured relative to atmospheric pressure 

Stagnation temperature t_ = 59° 

3 

The evaluation of the measured values was made according to the 
adiabatic 


T i ■ Mr 


K-l 


and the efflux equation 



with constant static pressure assured in the mixing region. 

Due to the imperfect readability of the thermoelement which, as 
mentioned before, is too rough for smaller temperature differences, it 
was impossible to measure the distribution of the stagnation tempera- 
tures over the mixing region. For the evaluation a linear drop of the 
stagnation temperatures along the mixing width was assumed. 

For the outer jet there results 

o _ 

t = 9 U A = ^1 me ^ ers P 61 * second 

for the jet issuing from the inner nozzle 

o — 

t^ = 13 Ui = 302 meters per second 

The inner jet therefore has, compared to the outer jet, an excess 

, O U 1 - Uf) 

temperature of 4 . For the velocity ratio the result = 0.5 

was obtained. 

Figure 17 shows the total pressure distribution 



$ 
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made dimensionless with the central value, for the various test sec- 
tions. The section near the nozzle which still shows the character of 
a turbulent pipe flow is represented in figure 18. Figure 19 shows, 
in addition, the variation of the total pressures along the jet axis. 

Figures 20 to 22 contain the corresponding representations for the 
velocities made dimensionless by the velocity u-^ of the jet issuing 
from the nozzle. 

■f 

As to the comparison with the theory, it must be noted that the 
velocity distribution at the exit from the nozzle is not rectangular, 
as assumed in the theory, but that it represents the profile of a 
turbulent pipe flow. (Compare fig. 21.) Hence it proves necessary to 
introduce the conception of the "effective diameter" in contrast to the 
geometric diameter. 

We define the effective nozzle diameter as the width of the 
rectangular velocity distribution of the amount u^ which is equiva- 
lent to the existing momentum distribution. That is, we calculate the 
effective nozzle diameter from the equation 

u l( u l - ^o) e g feCt ‘ = f Q "p - "o) T dr 

with the integral, which according to the theorem of conservation of 
momentum represents a constant, to be extended over an arbitrary cross 
section. 

In our case the integration over the cross section near the 
nozzle yields 


r effect . = 0 • 9^5r 


geom. 


Whereas the plotting over r\ - 


r - r 


geom. 


lets the test points appear 
r “ r effect. 


as still lying on one curve, the plotting over T] = — 

-A. 

results in a stagger of the velocity distributions with increasing — 

r 0 

toward negative T) . This stagger toward negative q expresses the 
immediately obvious fact that the isotacs of the flow field are curved 
toward negative tj (toward the jet axis) . 

x x k 

Figure 23 contains the theoretical curves for — =0 and — 

(the core end); in addition, the test points of the sections x = ?0 mil 
limeters and x = 45 millimeters were plotted. The agreement appears 
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to be good as far as the velocity gradient and the orientation in 
space in the central mixing region are concerned; the agreement in the 
transitions toward the jet core and the surrounding medium is less 
satisfactory. Deviations in these transitions are essential charac- 
teristics of the more recent Prandtl expression, hut are caused here 
probably mainly by the approximation character of our developments. 

For the core length there results according to the theory a value 
of x^ = 22.0r e f ; f ec ^, whereas the measurements along the jet axis 

(compare fig. 22 ) result in about x v = — ^9 = 20 . lr __ . . 

K 0.,9a 5 effect. 

It has to be noted that the experimental determination of the core end 
is affected by some uncertainty. 


IV. SUMMARY 


The spreading of a free rotationally symmetrical jet issuing from 
a nozzle represents a turbulent flow state. 

The theoretical investigation is based on the more recent Prandtl 
expression g = Kb J u^y - u m ^ n j for the momentum transport. The 

continuity equation and the equation of momentum are at disposal for 
calculation of the velocity distribution. In case of limitation to 

/ui - Uq 

small disturbances small quantity, where u-t is jet exit 

\ U 1 \ 

velocity, uq velocity of the surrounding medium J the equation of 
momentum may be linearized ' 


o . 

du / \/<3u Sul 

An approximate solution is constructed which is characterized by 
the fact that in boundary zones of the region as well as along the jet 
q = 0 in the interior of the region it has to be regarded as exact 
solution. 


/ 1 0 \ 

For arbitrary disturbances I arbitrary > 0 I the initial 

profile which corresponds to the velocity distributions of two mixing 
plane jets is determined by the fact that the transverse component in 
the jet core vanishes. The regularity found for small disturbances is 
taken as a basis for the further development of the profile from this 
initial profile.' 
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The transverse component of the flow is determined from the con- 
tinuity equation, with the use of the approximate function for the 
velocity component in the main flow direction. For simplicity a linear 
course of mixing width is assumed. 

The dimensions of the mixing region (limiting curve of the jet 
core d(x) and mixing width b(x)) are approximately calculated from 
the theorem of momentum 

p 00 

(u - Uq) vr - ^ / • u(u - U^jr dr (= rr X y) 

J r 

= rKb W ( u i - u o) ^ 

under assumption of a rectilinear course of the velocity distribution 

r - d 

over T), where q = — ^ — • 

In order to test the theory by experiment, a measurement was 

U-I - 

performed for — = 0.5 with a 5 millimeter nozzle. In order to 

U 1 

carry out the comparison with the theory, the conception of the effec- 
tive nozzle diameter is introduced which complies with the deviation 
of the effective velocity distribution for an issuing jet from the 
rectangular velocity distribution 

UiK - «o) = / o “ u(u - uo)r dr 

The agreement between theory and experiment is satisfactory. 



NACA TM 1311 


39 


V. REFERENCES 


1. Kuethe, Arnold M. : Investigations of the Turbulent Mixing Regions 

Formed by Jets. Jour. Appl. Mech., vol. 2, No. 3> 1935> 
pp. A87-A95. 

2. Prandtl, L.: Bemerkungen zur Theorie der freien Turbulenz. 

Z.f.A.M.M., Bd. 22, 1942. 

3. Gortler, H.: Berechnung von Aufgaben der freien Turbulenz auf Grund 

eines neuen Naherungsansatzes. Z.f.A.M.M., Bd.. 22, Nr. 5> 

Oct.. 1942, pp. 244-254. 

4. Tollmien, Walter: Berechnung turbulenter Ausbreitungsvorgange 

Z.f.A.M.M., Bd. 6 , Heft 6 , Dec. 1926, pp. 468-478. (Available 
as NACA TM IO 85 .) 

5. Pabst: Die Ausbreitungheisser Gasstrahlen in bewegter Luft. 

UM 8003, 1944. 

6 . Reichardt, H.: Uber eine neve Theorie der freien Turbulenz. 

Z.f.A.M.M., Bd. 21, 1941. 



appendix 


NO. 1. CALCULATION OF THE TRANSVERSE COMPONENT 



We substitute 


_u_ _ iN 

u i \T«\ u i 


pH 

/ 




d[] + \ 



□ 



u i - ^ 
u i 



du/ui _ J>0 
^2L = M U 1 

r O 



Su/ui ! /Uq 

""5n " ^\ u i 



e 
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If we assume a rectilinear course of the mixing width b(x), we 
have c-j_ * = 0 - 

We then obtain 



The evaluation of the integrals yields 
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where 


thus 


0 >) 






t _ 


_2_ 


J 0 


[] . 

e 


□ 


dG 


error integral 






G 

T) 


„ ^ U 1 - "0 
°1 T ' - °*3 6 ■ Ul - + 

Ui - llrv 

□ + 0.36 - L li r ^ - ct 2 


This results in 


j e -[f „ d , = n 1 °' 36j ^ ■ 02 


a i 


<1(0^) 


■ 7 ?f e ' n ‘ n dU 
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APPENDIX 


NO. 2. FOR CALCULATION OF THE DIMENSIONS OF THE CORE REGION 


The theorem of momentum with r = tq reads 


U - UQ ] / V 

\ U 1 / ' \ u l, 

X r=r 0 v 1 r =r Q 


a 

r 0 * SFI 


r drb d 7 

r Q "l\ U 1 / 




With the coordinate transfonnation T) 



we obtain 




r 0 



d u/u-^ 



- n(t>d' + b'd) 



dri 



u - 
U 1 


P^) 


(1 - Tj) 



We substitute 



h6 
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(bd' + b’d) 



dT} + dd 
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If one evaluates the integrals, one obtains 


' a - V l/y \ m / U 1 - ”0 

v U 1 Mr.T 0 T ° ' \ U 1 





= bb 1 


f u i - " o 
V u i 



ir(l - Tl)dTJ + 


r cr d - 


^o h, - 

u i\ u i 


2 , 

t) drj 


r 0 -d 


+ (b’d + bd') 


^i^Uq 
Ui 



nO- 


r 0 -d 


V U 1 • Vi r 1 


u i\ u i 


T1 *1 


r 0 -d 


+ dd» 


/u i ~ u o\ r 1 


Ui 


r 0 -d 


(1 - T] ) dt) + 



T))dT1 + 


V } 2 1 

— j 5- 





k8 


NACA TM 1311 


If we insert these expressions into the equation of momentum and 
order, we obtain 


. bb' 


'ifSr'f M’ 1 


/ Ut - 


"0 


U 1 


iN’fM ■ * 


1 Up/ U i - V ,/ r 0 ~ d \ 3 
3 U X \ U X /\ b ) 


+ (bd’ + b’d) 


l/ u l - ^of/ro - d 

2 l U 1 ) \ * 






2 / U 1 - « 0\ 2 / r 0 * d \ 3 1 V *! ' M . 1 ^o/ u l - u p \/ r Q - d 

n - ^ - ) r ^ / " 5 2 v u i A b 
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which finally leads to 



Translated by Mary L. Mahler 
National Advisory Committee 
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Figure 2.- Auxiliary quantities for calculation of the velocity 

distribution. 






Figure 3. - Transverse component distribution. 
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Figure 4.- Transverse component distribution. 
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Figure 5." Transverse component distribution. 
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Figure 6.- Transverse component distribution. 
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Figure 8.- Dimensions of the core region. 
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Figure 9.- Gore lengths. 
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Figure 11.- Mixing width variation. 
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Figure 13.- Parameter value of the plane jet rim. 
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Figure 14.- k as a function of c 
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Figure 15.- <tq as a function of c. 




66 


NACA TM 1311 























c 






+ 









E 

a 






X 

rv 









• 

i5 






+ 




i? 



* 
__ H 

qo Si 2 

4 + 0 

CO 

II 




X 




— s 
■ 

-J 





0 




~8 

+ 















X 















3 

+ 

















X 















♦ 

ri 

















1 

4 













□ 

O 

♦ 

CM- 

X 



— F' — 

s 









H 

♦x” 1 

1 * 




— n 

m 






X 

_rt_ 

X 



« 







“1 




■' j 

< 

+ 

—S 

0 

0 

u 1 




l 






4 

f + 

3 + 














1 

.4 

IT — 

f 0 

0 



■ -~l 










| 











|| 



I 











c 














“ ^ 
1 

3 

.s 



□ 

( 














X 


8^ 

+ 
















K 

O 

< + 


'l— 















X 























I 

-c 














l*3i? t 
+ + ; 
i4*i£ 

j 1 

5 




w 

r> 







■ - «P 






C 

3 

J 







co 


6 

6 




o 


*0 


* 


■*> 


<NJ 


O 


04 

t 


n 


/ 


Figure 17. 
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Figure 18. 



68 


NACA m 1311 



Figure 19 
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Figure 20 
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Figure 23. - Comparison between theory and experiment. 




